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Abstract. - In order to characterize populations of Diplodus vulgaris (Geoffroy St.-Hilaire, 1817) along the 
Tunisian coasts and to test the prospective role of geographical barrier of the Siculo-Tunisian Strait, we per¬ 
formed a study of morphological and genetic variability of samples collected in six sites distributed on each side 
of the Strait. Two types of morphological descriptors, 25 truss elements and six complementary conventional 
linear measurements, and genetic markers (allozymes of four enzyme systems) were used. Multivariate analyses 
revealed morphological variations among the prospected localities. El Biban sample differs by the suborbital 
part of the head and the eye diameter, while those of Ghar El Melh and Mahdia diverge mainly by the length and 
height of the head. Compared to all remaining samples, El Biban specimens showed a large head. The study of 
genetic polymorphism showed the absence of any genetic structuring among Diplodus vulgaris and the presence 
of an homogeneous (global FST = 0.004) panmictic unit (f = 0.046 NS). Thus, the Siculo-Tunisian Strait does 
not seem to act as a barrier between the populations, at least at the scale of the Tunisian coast. The morphological 
diversity observed seems to reflect the environmental heterogeneity and is not structured by the Siculo-Tunisian 
Strait (phenotypic plasticity). 
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Resume. - Variations morphologique et genetique de Diplodus vulgaris le long des cotes tunisiennes. 

Pour caracteriser les populations de Diplodus vulgaris (Geoffroy St.-Hilaire, 1817) le long des cotes tuni¬ 
siennes et tester le role eventuel de barriere geographique du detroit siculo-tunisien, nous avons etudie la variabi- 
lite morphologique et genetique d’echantillons collectes dans six sites de part et d’autre du detroit. Deux types de 
descripteurs morphologiques, 25 elements de truss et six mesures lineaires conventionnelles complementaires, 
et des marqueurs genetiques (allozymes de quatre systemes enzymatiques) ont ete utilises. Des analyses multi- 
variees ont revele des variations morphologiques entre les localites prospectees. Ainsi, l’echantillon d’El Biban 
differe par la partie suborbitale de la tete et le diametre de l’ceil alors que ceux de Ghar El Melh et de Mahdia 
divergent principalement par la longueur et la hauteur de la tete. Par rapport a Tensemble des echantillons, les 
specimens d’El Biban ont une plus grosse tete. L’etude du polymorphisme genetique a montre l’absence de toute 
structuration genetique chez Diplodus vulgaris et la presence d’une unite panmictique (f= 0.046 NS) homogene 
(f'STglobal = 0-004). Ainsi, le detroit siculo-tunisien ne semble pas jouer le role de barriere entre les populations, 
au rnoins a l’echelle de la cote tunisienne. La diversite morphologique (plasticite phenotypique) semble refleter 
Theterogeneite environnementale et n’est pas structuree par le detroit siculo-tunisien. 


Marine habitats are dispersive systems, which do not 
exhibit clear physical barriers, favouring gene exchange 
among populations (Palumbi, 1994) and contributing to 
limit the action of genetic differentiation factors among 
populations (Turan et al., 2004). However, some geograph¬ 
ic marine areas are recognized to play the role of transition 
zone for several marine species (Borsa et al., 1997; Bargel- 
loni et al., 2005; Patarnello et al., 2007; Galarza et al., 2009; 
Mejri et al., 2012). In the Mediterranean Sea, the Siculo- 
Tunisian Strait (STS), located between Cap Bon in Tunisia 
and Mazara del Vallo in Italy, constitutes a physical bound¬ 
ary between the eastern and western Mediterranean basins 
(Beranger et al., 2004). It seems to constitute a geographic 
barrier for some marine fish species, which are subject to 
genetic transition in this area (Bahri-Sfar et al., 2000; Ste- 


fanni and Thorley, 2003; Zardoya et al., 2004; Mejri et al., 
2009, 2011). Several genetic studies, especially on sparids, 
were carried out along this area, highlighting different genet¬ 
ic patterns (Ben Slimen et al., 2004; Hammami et al., 2007; 
Fassatoui et al., 2009; Kaoueche et al., 2011). The transition 
between the two basins seems to be perceived in different 
ways depending on the species since a strong to moderate 
genetic divergence is observed in some species, whereas 
others show homogenous pattern. For instance, for Lithog- 
nathus mormyrus , no evidence of genetic structuring was 
observed (Hammami et al., 2007). In contrast, a strong dif¬ 
ferentiation was detected among Tunisian Sparus aurata 
between groups located in the eastern and western Mediter¬ 
ranean basins (Ben Slimen et al., 2004), whereas a slight 
differentiation was observed between distant samples of 
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Pagellus erythrinus (Fassatoui et al., 2009). Finally, in the 
case of Diplodus sargus , a differentiation was revealed by 
allozymes and was attributed to small-scale, not patterned, 
genetic heterogeneity among local populations, whereas the 
mtDNA data indicate a genetic homogeneity of the white 
seabream populations among the same localities across the 
STS (Kaoueche et al., 2011). 

In addition, a sufficient degree of isolation may result in 
notable morphological differentiation among stocks within 
a species (Pepin and Carr, 1993). Beside the primary inter¬ 
est in understanding the patterning of fish populations, the 
description of the morphological variation between popula¬ 
tions and the identification of correlated degree of genetic 
pattern help to implement the appropriate management 
strategies of Mediterranean fisheries. In this context, mor¬ 
phometric characters (Turan and Basusta, 2001; Palma and 
Andrade, 2002; Turan, 2004; Bahri-Sfar and Ben Hassine, 
2009; Hammami et al., 2011) or morphometric characters 
associated to genetic data (Cabral et al., 2003; Mejri et al., 
2012) have been commonly used as markers to characterize 
and delimit geographic stocks (Arculeo et al., 2003). 

The two-banded seabream, Diplodus vulgaris is one 
of the most frequent and abundant sparid in the Mediterra¬ 
nean Sea and along the eastern Atlantic coast (Bauchot and 
Hureau, 1986). The species is generally found associated 
with both rocky and sandy bottoms. luveniles inhabit coastal 
lagoons and estuaries (Macpherson, 1998). D. vulgaris life 
cycle is characterized by an offshore pelagic larval disper¬ 
sion during 29-58 days followed by a nectobenthic adult 
phase. These characteristics give great dispersal potential, 
since in marine environments, early-life-history traits and 
dispersal capabilities may influence exchanges of gene flow 
between populations (Palumbi, 2004; Galarza et al., 2009). 

In this context, we have undertaken the characteriza¬ 
tion of Mediterranean populations of two-banded seabream 
across the STS. Effectively, despite the commercial impor¬ 
tance of this species, studies are scarce. In this regard, a first 
genetic study of Mediterranean Diplodus vulgaris popula¬ 
tions, using allozyme markers in seven samples from Tyrrhe¬ 
nian, Aegean, Ionian and Adriatic Seas and also both sides 
of the STS, demonstrated a distinction of the Adriatic Sea 
sample only (Arculeo et al., 2003). However, Galarza et al. 
(2009) have shown a genetic transition across the Almeria- 
Oran front using microsatellite markers, but with sampling 
that did not cover the entire distribution area of D. vulgaris 
in the Mediterranean Sea. In the present study, we aimed to 
check whether morphological variation in D. vulgaris was 
in relation with the STS, then we examined the possible cor¬ 
relation with genetics and ecological features of the studied 
populations. For that purpose, we collected samples along 
the Tunisian coasts, on both sides of the STS in order to 
characterize the population of D. vulgaris in this transition 
area with heterogeneous environments. The morphometric 
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Figure 1. - Locations of sampling sites along the Tunisian coasts. 
★ : marine samples. •: lagoon samples. STS: Siculo-Tunisian 
Strait. West-Med: Western Mediterranean basin. East-Med: Eastern 
Mediterranean basin. 

study was based on body shape descriptors (conventional 
linear and truss measurements), while genetic investigation 
was based on allozyme polymorphism. 

MATERIAL & METHODS 

Sampling 

Samples were collected from six study sites (4 marine 
and 2 lagoonal), distributed on each side of the STS and cov¬ 
ering northeastern (Bizerta Bay, Ghar El Melh lagoon and 
Goulette) and southeastern (Mahdia, Gabes Gulf, and El 
Biban Lagoon) sectors of Tunisian coastline (Fig. 1). These 
localities differ by their environmental features (Tab. I). 
Sample size ranged from 27 to 51 individuals (Tab. I). All 
fishes were adults. They were captured using trammel nets 
during the period between September 2007 and December 
2008. The samples were brought to the laboratory in order 
to take the defined morphometric measurements, dissect and 
remove pieces of muscle and liver, which we stored at -20°C 
for allozyme exams. 

Morphometric data 

In order to describe the shape of collected fish, the truss 
approach was used to create a network describing fish body 
outline (Strauss and Bookstein, 1982). This approach is 
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Table I. - Sample locations of specimens of D. vulgaris and main environmental features of study sites. MSL: mean standard length. 
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Figure 2. - A: Location of the 11 landmarks (1-11) used for con¬ 
structing the truss network on D. vulgaris and the six additional 
points (12-17) used to draw the conventional linear measurements. 
Landmarks and additional points illustrated as black dots. Truss 
network illustrated as continuous lines. Conventional linear meas¬ 
urements illustrated as discontinuous lines. B: Discriminated head 
region and discrimination related variables. 

optimal for extracting body shape differences among popu¬ 
lations. The left side of each fish was photographed with a 
high-resolution digital camera mounted on a tripod. Mor¬ 
phological landmarks were selected around the outline of the 
fish form (Fig. 2). Landmark coordinates were performed 
on digital images using Visilog 6.480 software. Eleven 
landmarks determining 25 measurements were recorded 
(Fig. 2). Six points were added to obtain complementary 
conventional linear measurements, i.e. eye diameter, upper 
and lower head height, snout length, operculum length, and 
width of the pectoral fin base (Fig. 2; Annex I). Thus, 31 
measurements were analysed (Tab. I). In order to eliminate 
any size effect, size-dependent variation of morphometric 
measurements was removed using an allometric transforma¬ 
tion (Reist, 1985): M tran s = log M- (3 (log SL- log SL mean ) 
where M trans is the transformed measurement, M the original 
measurement, [3 the within-group slope regressions of log M 
against log SL, SL the standard length of the fish and SL mean 
the overall mean of the standard length. 

The analysis of variance (ANOVA) was performed to test 
the statistical significance of mean variations among samples 
and for each measurement. Multivariate analyses were used 
to elucidate pattern of the morphological variation between 
samples. The t-test was established to infer whether the aver¬ 
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ages of each variable are significantly different between two 
considered samples. 

The discriminant function analysis (DFA) was assessed 
in order to carry out the degree of similarity between the 
studied samples and the relative importance of each meas¬ 
urement for group separation. Wilks’ X values were esti¬ 
mated to test the significance of observed discrimination for 
a combination of variables. The classification success rate 
(PCS) was evaluated based on the percentage of individuals 
correctly assigned into the original sample. Statistical analy¬ 
ses were performed using R 2.11.1 software. 

Allozyme data 

Protein extraction was released in Tris buffer (pH = 6.8). 
Electrophoresis and specific staining were performed 
according to Pasteur et al. (1987). Four enzymatic systems, 
showing clear zymograms, were used in this analysis: alco¬ 
hol deshydrogenase EC 1.1.1.1 (ADH; TC 8), isocitrate 
deshydrogenase EC 1.1.1.4.2 (ICD; TC 8), lactate deshydro¬ 
genase EC 1.1.1.27 (LDH; TC 8), glucose phosphomutase 
EC 5.4.2.1 (PGM; TC 8) and glucose phosphate isomerase 
EC 5.3.1.9 (PGI; TC8). Eight polymorphic loci were scored 
(ADH*, ICD-1*, ICD-2*, LDH*, PGI-1*, PGI-2*, PGI-3* and 
PGM*). Nomenclature of loci and alleles was written accord¬ 
ing to Shaklee et al. (1990). 

Genetic diversity was described through estimating the 
unbiased expected heterozygosities ( H e ) (Nei and Chesser, 

1983) , mean number of alleles per locus (A m ), and number 
of exclusive alleles (A e ). Departure from Hardy-Weinberg 
Equilibrium (HWE) was tested with the estimator/of Wright 
fixation index (Eis) (Wright, 1969) according to Weir and 
Cockerham (1984). The estimator 9 (Weir and Cockerham, 

1984) of Fst was calculated between populations in order 
to estimate the genetic differentiation between samples. All 
analyses were carried out using Genetix v.4.05 software 
(Belkhir et al., 2004). 


RESULTS 

Morphometric analysis 

The univariate analysis (ANOVA) shows significant 
differences (P < 0.001) among localities for all studied 
variables (Tab. I). The projection of the six samples on 
the plan defined by DF1 and DF2 illustrates the variation 
among samples (Fig. 3). More precisely: the first function 
DF1, which explains 28% of the global variation, high¬ 
lights the distinction of the southern lagoon sample (El 
Biban lagoon), whereas the second function (DF2), which 
explains 24% of the global variation, allows the discrimina¬ 
tion of the northern lagoon (Ghar El Melh lagoon) and an 
eastern marine sample, namely Mahdia; the remaining three 
localities are undifferentiated according to DF1 and DF2 
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(Fig. 3). The discrimination of El Biban lagoon is essen¬ 
tially highlighted by two truss measurements located in 
the head region and the eye diameter (VI, V2, V30; Fig. 2; 
Annex I). In fact, the head of El Biban lagoon specimens has 
larger eyes and a longer predorsal edge than all other sam¬ 
ples. The application of t-test among these three discrimi¬ 
native measurements showed significant values (Tab. II). 
El Biban lagoon sample has the highest means of the truss 
variables (V2 B ;t, an = 2.261 mm; V25 B ;b an = 6.984 mm) and 
the eye diameter (V30 B ;b an = 1 -274 mm) than Ghar El Melh 
lagoon (V2ch a rM= 1-319 mm, V25charM = 4.697 mm and 
V30 G h ar M = 0.828 mm) and Mahdia (V2 Mahdia= 1.054 mm, 
V25 Ma hdia = 4.837 mm and V30 Ma hdia= 1-100 mm). 

Ghar El Melh lagoon sample and Mahdia marine sample 
are discriminated along the second function DF2, defined by 
three measurements related to the head dimensions (Fig. 2: 
VI, V12, V29; Tab. I). Significant Average differences 
between Ghar El Melh and Mahdia are highlighted by V1 
and V29 (Tab. II) with the lowest average for the front head 
length (Fig. 2: VI) in Ghar El Melh population. In addition, 
pairwise significant average comparisons show that in all 
cases, El Biban lagoon fish have the biggest head (Tab. II). 


Figure 3. - DFA scores of morphomet¬ 
ric characters using conventional linear 
measurements and truss elements on 
the plan DF1-DF2. 

Observed morphometrical distinction among Tunisian D. 
vulgaris population was confirmed by a significant value of 
Wilks' X criterion (Wilks X = 0.0236; F = 7.084; P < 0.001). 

Finally, the overall assignment of individuals into their 
original sample (PCS) was 82.16%. The percentage values of 
specimens correctly classified into the original groups were 
ranged from 61% to 100% (Tab. Ill), the highest percentage 
being observed in El Biban and Mahdia samples. This classi¬ 
fication confirmed the clear discrimination of these samples 
(Fig. 3). 

Allozymes 

Twenty-one alleles were observed among the eight 
typed polymorphic loci. Four unique alleles were detected 
(Tab. IV). Three of them are observed in the southern coasts 
(ADH* 120 in El Biban lagoon; PGI-2* 120 and PGI-3* 120 
in Gabes Gulf). A fourth one is detected in Ghar El Melh 
lagoon (PGI-3* 80). However, the all four are considered as 
rare alleles. Allele frequencies ranged from 0.019 for PGI-3* 
80 to 0.027 lav ADH* 120. Diversity parameters showed that 
the highest values of average number of alleles (A m = 2.125) 
and average expected heterozygosity (He = 0.072 ± 0.014) 
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are observed in the sample from the Gabes Gulf (Tab. IV). 

Global Fis value estimated for all samples is low and not 
significant (f= 0.046). Multilocus/values in each locality 
are not significant (Tab. IV). They range from -0.024 NS in El 
Biban lagoon sample to 0.209 NS in Goulette sample. 

Global Fsr value is also low and not significant (0.004), 
showing the absence of any differentiation among studied 

Table II. - Average comparisons using Student’s t-test; ddl: degree 
of liberty. A: t-test values calculated between discriminated vari¬ 
able averages for distinctness samples along DF1. B: t-test values 
calculated between discriminated variable averages for distinctness 
samples along DF2. See table I for sample site codes. 


A 

ddl 

V2 

V25 

V30 

LBIBV/MBIZV 

84 

16.118*** 

9.098*** 

45.728*** 

LBIBV/LGEMV 

63 

65.527*** 

43.627*** 

141 47 *** 

LBIBV/MGOUV 

69 

15.374*** 

8.452*** 

87.418*** 

LBIBV/MMAHV 

71 

92.277*** 

26.923*** 

29.823*** 

LBIBV/MGABV 

87 

31.302*** 

19.283*** 

76.911*** 


B 


Student’s t-test 


ddl 

VI 

V12 

V29 

MMAHV/MBIZV 

82 

-4.595*** 

-6.226*** 

-16.551*** 

MMAHV/LGEMV 

61 

-4.566* 

ns 

-7.817* 

MMAHV/MGOUV 

67 

ns 

ns 

-14.329** 

MMAHV/MGABV 

85 

8.174*** 

4.726** 

9.115* 

MMAHV/LBIBV 

71 

-14.577*** 

-18.216*** 

-58.691*** 

LGEMV/MBIZV 

74 

-9.571* 

-10.129*** 

-21.524*** 

LGEMV/MGOUV 

59 

ns 

-4,994** 

-21.014*** 

LGEMV/MGABV 

77 

ns 

ns 

ns 

LGEMV/LBIBV 

63 

-23.68*** 

-34.769*** 

- 68 . 666 *** 
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samples. The estimation of pairwise 0 values confirms this 
homogeneity. Pairwise 0 values are not significant and range 
from 0.000 to 0.013 between the Goulette sample, located 
in the north Tunisian coast, and the El Biban lagoon sample 
which is a southern sample (Tab. V, Fig. 1). 

The estimation of pairwise genetic distances between 
the samples reveal low values (Tab. V), improving obtained 
results with 0 estimator. 


DISCUSSION 

The analyses of morphometric data of Diplodus vulgaris 
collected in marine and lagoonal environments along the 
Tunisian coasts on both sides of the STS revealed morpho¬ 
logical variation among localities but do not show a relation 
with the presence of the STS. The lagoon sample at Ghar 
El Melh and at El Biban and the southeast marine sample 
at Mahdia exhibit each distinguishing head morphology, 
whereas all other marine samples exhibit close median mor¬ 
phology of the head (Bizerta, Goulette and Gabes Gulf). The 
specimens from El Biban lagoon are notably characterized 
by their biggest head when compared to other samples. 

Morphological differences in the head region have been 
documented in several marine fish species (Turan, 2004; 
Wainwright et al., 2004; Mejri et al., 2012) and similar 
results were obtained in other sparid species. For instance, a 
significant degree of morphological dissimilarities was evi¬ 
denced between sparid populations of Diplodus sargus, D. 
puntazzo from four Mediterranean (Spain, Italy and Greece) 
and Atlantic (Portugal) localities, among which more than 


Table III. - Percentage of specimens classified in each group. 



El Biban 
lagoon 

Ghar El Melh 
lagoon 

Bizerte 

Bay 

Gabes 

Gulf 

Goulette 

Mahdia 

El Biban lagoon 

100 

- 

- 

- 

- 

- 

Ghar El Melh lagoon 

- 

82 

4 

7 

7 

- 

Bizerte Bay 

10 

2 

70 

6 

10 

2 

Gabes Gulf 

6 

- 

14 

80 

- 

- 

Goulette 

3 

- 

34 

3 

61 

- 

Mahdia 

- 

- 

- 

- 

- 

100 


Table IV. - Genetic diversity parameters in D. vulgaris samples from the Tunisian coasts: sample size (AO, number of exclusive alleles (A e ), 
mean number of alleles per locus and per population (A m ), average unbiased expected heterozygosity (7/ e ), Weir and Cockerham (1984) 
multilocus fixation index (/) measuring departure from theoretical Hardy-Weinberg equilibrium, confidence interval (±), non significant 
(NS). 



Bizerte Bay 

Goulette 

Ghar El Melh 
lagoon 

Mahdia 

Gabes Gulf 

El Biban lagoon 

N 

48 

33 

27 

35 

51 

37 

Ae 

- 

- 

1 

- 

2 

1 

Am 

2 

1.625 

2 

1.5 

2.125 

1.625 

He 

0.062 ±0.014 

0.056 ±0.022 

0.066 ±0.016 

0.030 ±0.013 

0.072 ±0.014 

0.035 ±0.013 

f 

-0.047NS 

0.209 NS 

0.112 NS 

-0.030 NS 

0.091 NS 

-0.024 NS 


116 


Cyhium 2013, 37(1-2) 











































Kaoueche etal. 


Morphology and genetic o/'Diplodus vulgaris on Tunisian coasts 


Table V. - Genetic distance values (Nei, 1987) (above diagonal) and 9 values (below diagonal), estimator of Fst according to Weir and 
Cockerham (1984) among D. vulgaris samples for allozymic data. 



Mahdia 

Bizerte Bay 

Goulette 

Ghar El Melh 
lagoon 

El Biban lagoon 

Gabes Gulf 

F st 

0 

0 

0 

0 

0 

0 

Mahdia 


0.000 

- 0.000 

0.000 

0.000 

0.000 

Bizerte Bay 

0.00 INS 


0.000 

0.000 

0.000 

0.001 

Goulette 

0.000 NS 

0.000 NS 


0.000 

0.001 

0.001 

Ghar El Melh lagoon 

0.011 NS 

0.006 NS 

0.002 NS 


- 0.000 

0.000 

El Biban lagoon 

0.006 NS 

0.002 NS 

0.013 NS 

0.003 NS 


0.000 

Gabes Gulf 

0.005 NS 

0.008 NS 

0.010 NS 

0.001 NS 

0.007 NS 



50% were head-related (Palma and Andrade, 2002). More¬ 
over, a morphological discrimination, assigned to the head 
length allows the distinction between lagoon populations 
of Lithognathus mormyrus (Hammami et al., 2011,2013). 
Finally, similar results regarding the head morphology were 
obtained by Sara et al. (1999) on cultivated D. puntazzo 
reared under different conditions. So, if well documented, 
in the Mediterranean region, such morphological discrep¬ 
ancies in fish populations are attributed either to the impact 
of habitat differences (Palma and Andrade, 2002; Cabral et 
al., 2003; Turan, 2004; Hammami et al., 2011; Kaoueche, 
2012), or to the interaction of both genetic and environmen¬ 
tal components (Favaloro and Mazzola, 2006; Bahri-Sfar 
and Ben Hassine, 2009; Mejri et al., 2012), with the rupture 
or reduction of gene flow attributed to the existence of fronts 
and straits acting as barriers (Bahri-Sfar et al., 2000; Bargel- 
loni et al., 2005; Galarza et al., 2009; Bahri-Sfar and Ben 
Hassine, 2009; Zittari-Chatti et al., 2009; Mejri et al., 2012). 
In this study, genetic analyses based on allozyme markers 
demonstrates the homogeneity among our samples. There¬ 
fore, we can consider that Tunisian Diplodus vulgaris con¬ 
stitute a homogeneous panmictic unit and that STS does not 
seem to act as a barrier to gene flow for this species at the 
scale of the Tunisian coast. This can be explained by a quite 
long pelagic larval duration of D. vulgaris and larval disper¬ 
sal (29-58 days; Bauchot and Hureau, 1986), which favour 
gene flow and population homogenization. These results are 
not so surprising since low level of genetic structure was 
highlighted in the Mediterranean region for D. vulgaris, with 
only differentiation across the Almeria-Oran Front (Galar¬ 
za et al., 2009) and the discrimination of the Adriatic Sea 
population (Arculeo et al., 2003). Moreover, a similar result 
(homogeneity) was observed with the related species D. sar- 
gus (Kaoueche et al., 2011), which shares life history traits 
but displays a shorter phase of larval dispersal (21-30 days). 

So, the morphological variations observed between the 
Tunisian two-banded seabream cannot be attributed to genet¬ 
ic differentiation. It is thus related to environmental differ¬ 
ences between the studied sites and could reflect an adaptive 
divergence of populations. 


The fact that the morphological variations concentrate in 
the head region leads us to search for eventual causes in fish 
diet (e.g. Walls et al., 1993; Palma and Andrade, 2002; Ham¬ 
mami et al., 2013). Indeed, Diplodus vulgaris is a demersal 
fish that feeds on benthic invertebrates (Fischer et al., 1987). 
However, depending on environmental characteristics pre¬ 
vailing in each of these sites, the availability and abundance 
of potential prey vary from one site to another (Rosecchi, 
1985). On the east coast of Tunisia, their preferential prey 
are crustaceans, molluscs, annelids and echinoderms, while 
secondary preys correspond to sponges and platyhelminths 
(Tekaya, 2009). In contrast, in the Gulf of Gabes, D. vulgaris 
feeds preferentially on bivalves and crustaceans, and sec¬ 
ondarily on annelids and cephalopods (Bradai, 2000). This 
difference in diet composition may reflect the prey availabil¬ 
ity in relation to the heterogeneity of target environments. 
However, direct correlation between head shape of a popula¬ 
tion and its environment cannot be assumed. This is prob¬ 
ably because of mainly multifactorial controls on prey avail¬ 
ability. Indeed, precise information on the diet of the given 
populations (stomach contents along several seasons) is nec¬ 
essary to explain further head morphology. 

So, since it is admitted that when a species colonizes a 
range of heterogeneous environments, a single phenotype 
is unlikely to be associated with high fitness throughout the 
range (Via et al., 1995), a species must adapt to succeed, 
notably by being able to match phenotypes with local envi¬ 
ronments (J0rgensen et al., 2008). Such adaptations may 
imply changes in morphology, reproduction or survival fea¬ 
tures (Stearns, 1983). Our study suggests that, on Tunisian 
coasts, the two-banded seabream is able to adapt to local 
environments through phenotypic plasticity alone, which is 
considered an optimal strategy to coping in heterogeneous 
environments (J0rgensen et al., 2008). Conversely, the STS 
does not play any role in the patterning of the populations of 
Diplodus vulgaris, at least at the scale of the Tunisian coast. 
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Annex I. - Loadings from discriminant function analysis. Univariate analysis (ANOVA) for all measurements. Significance levels: 
*: P < 0.05; **: P < 0.01; ***: P < 0.001. %: discrimination percentage of each function. In bold: more discriminant variables along the 
two axes. 
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0.363 

V26: Xl_12 

17.997*“ 

0.011 

0.475 

V27: Xl_15 

16.636*“ 

-0.493 

0.272 

V28: X2_15 

17.316*“ 

-0.371 

0.484 

V29: 11_15 

27.932*“ 

-0.357 

0.506 

V30: X13_14 

24.893*** 

-0.647 

0.242 

V31: X16 17 

22.901*** 

-0.542 

0.276 
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